{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}C:\Stata\log_figure1.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res}14 Jul 2018, 19:50:02
{txt}
{com}. 
. set more off
{txt}
{com}. set scheme s1manual
{txt}
{com}. 
. probit civicnrt1 i.smlmegip##i.highdemoc unrest1 majpow1 lnCapRatio alliance i.highdemoc##i.highdemoc2 tradedep NonCNRyrst1 NonCNRyrs2t1 NonCNRyrs3t1 CNRyrst1 CNRyrs2t1 CNRyrs3t1, cluster(dyadid), if Contig==1

{res}{txt}Iteration 0:{space 3}log pseudolikelihood = {res: -4880.502}  
Iteration 1:{space 3}log pseudolikelihood = {res:-794.92428}  
Iteration 2:{space 3}log pseudolikelihood = {res:-732.06632}  
Iteration 3:{space 3}log pseudolikelihood = {res:-663.41764}  
Iteration 4:{space 3}log pseudolikelihood = {res:-660.65943}  
Iteration 5:{space 3}log pseudolikelihood = {res:-660.52318}  
Iteration 6:{space 3}log pseudolikelihood = {res:-660.52277}  
Iteration 7:{space 3}log pseudolikelihood = {res:-660.52277}  
{res}
{txt}Probit regression{col 49}Number of obs{col 67}= {res}    34,962
{txt}{col 49}Wald chi2({res}16{txt}){col 67}= {res}   1982.36
{txt}{col 49}Prob > chi2{col 67}= {res}    0.0000
{txt}Log pseudolikelihood = {res}-660.52277{txt}{col 49}Pseudo R2{col 67}= {res}    0.8647

{txt}{ralign 86:(Std. Err. adjusted for {res:967} clusters in dyadid)}
{hline 21}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 22}{c |}{col 34}    Robust
{col 1}           civicnrt1{col 22}{c |}      Coef.{col 34}   Std. Err.{col 46}      z{col 54}   P>|z|{col 62}     [95% Con{col 75}f. Interval]
{hline 21}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 10}1.smlmegip {c |}{col 22}{res}{space 2} .0215505{col 34}{space 2} .1076076{col 45}{space 1}    0.20{col 54}{space 3}0.841{col 62}{space 4}-.1893564{col 75}{space 3} .2324575
{txt}{space 9}1.highdemoc {c |}{col 22}{res}{space 2}-.0601124{col 34}{space 2} .1578793{col 45}{space 1}   -0.38{col 54}{space 3}0.703{col 62}{space 4}-.3695502{col 75}{space 3} .2493254
{txt}{space 20} {c |}
{space 2}smlmegip#highdemoc {c |}
{space 16}1 1  {c |}{col 22}{res}{space 2} .4561446{col 34}{space 2}  .206608{col 45}{space 1}    2.21{col 54}{space 3}0.027{col 62}{space 4} .0512005{col 75}{space 3} .8610888
{txt}{space 20} {c |}
{space 13}unrest1 {c |}{col 22}{res}{space 2}  .002265{col 34}{space 2} .0079695{col 45}{space 1}    0.28{col 54}{space 3}0.776{col 62}{space 4} -.013355{col 75}{space 3}  .017885
{txt}{space 13}majpow1 {c |}{col 22}{res}{space 2} .1484077{col 34}{space 2} .1799964{col 45}{space 1}    0.82{col 54}{space 3}0.410{col 62}{space 4}-.2043788{col 75}{space 3} .5011942
{txt}{space 10}lnCapRatio {c |}{col 22}{res}{space 2} .1007282{col 34}{space 2} .0430348{col 45}{space 1}    2.34{col 54}{space 3}0.019{col 62}{space 4} .0163814{col 75}{space 3} .1850749
{txt}{space 12}alliance {c |}{col 22}{res}{space 2} .0871493{col 34}{space 2} .0908383{col 45}{space 1}    0.96{col 54}{space 3}0.337{col 62}{space 4}-.0908906{col 75}{space 3} .2651892
{txt}{space 8}1.highdemoc2 {c |}{col 22}{res}{space 2} .0502196{col 34}{space 2} .1355085{col 45}{space 1}    0.37{col 54}{space 3}0.711{col 62}{space 4}-.2153721{col 75}{space 3} .3158114
{txt}{space 20} {c |}
highdemoc#highdemoc2 {c |}
{space 16}1 1  {c |}{col 22}{res}{space 2}-.4327521{col 34}{space 2} .2316263{col 45}{space 1}   -1.87{col 54}{space 3}0.062{col 62}{space 4}-.8867313{col 75}{space 3} .0212271
{txt}{space 20} {c |}
{space 12}tradedep {c |}{col 22}{res}{space 2}  .060609{col 34}{space 2} .0225428{col 45}{space 1}    2.69{col 54}{space 3}0.007{col 62}{space 4} .0164259{col 75}{space 3} .1047921
{txt}{space 9}NonCNRyrst1 {c |}{col 22}{res}{space 2}-.3194074{col 34}{space 2} .0323339{col 45}{space 1}   -9.88{col 54}{space 3}0.000{col 62}{space 4}-.3827806{col 75}{space 3}-.2560342
{txt}{space 8}NonCNRyrs2t1 {c |}{col 22}{res}{space 2} .0131616{col 34}{space 2} .0020979{col 45}{space 1}    6.27{col 54}{space 3}0.000{col 62}{space 4} .0090498{col 75}{space 3} .0172733
{txt}{space 8}NonCNRyrs3t1 {c |}{col 22}{res}{space 2}-.0001578{col 34}{space 2} .0000351{col 45}{space 1}   -4.50{col 54}{space 3}0.000{col 62}{space 4}-.0002265{col 75}{space 3} -.000089
{txt}{space 12}CNRyrst1 {c |}{col 22}{res}{space 2} .5455503{col 34}{space 2} .0399733{col 45}{space 1}   13.65{col 54}{space 3}0.000{col 62}{space 4} .4672041{col 75}{space 3} .6238966
{txt}{space 11}CNRyrs2t1 {c |}{col 22}{res}{space 2}-.0247239{col 34}{space 2} .0031303{col 45}{space 1}   -7.90{col 54}{space 3}0.000{col 62}{space 4} -.030859{col 75}{space 3}-.0185887
{txt}{space 11}CNRyrs3t1 {c |}{col 22}{res}{space 2} .0003119{col 34}{space 2} .0000575{col 45}{space 1}    5.43{col 54}{space 3}0.000{col 62}{space 4} .0001993{col 75}{space 3} .0004245
{txt}{space 15}_cons {c |}{col 22}{res}{space 2}-1.007781{col 34}{space 2}  .146503{col 45}{space 1}   -6.88{col 54}{space 3}0.000{col 62}{space 4}-1.294921{col 75}{space 3}-.7206401
{txt}{hline 21}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{p 0 6 2}Note: 1306 failures and 22 successes completely determined.{p_end}

{com}. 
. replace smlmegip=1
{txt}(35,233 real changes made)

{com}. replace highdemoc=1
{txt}(30,095 real changes made)

{com}. replace unrest1=1
{txt}(42,211 real changes made)

{com}. replace majpow1=0
{txt}(3,963 real changes made)

{com}. replace lnCapRatio=-1.936979
{txt}(47,606 real changes made)

{com}. replace alliance=0
{txt}(21,990 real changes made)

{com}. replace highdemoc2=0
{txt}(22,049 real changes made)

{com}. replace tradedep=.2491212
{txt}(47,606 real changes made)

{com}. replace NonCNRyrst1=0
{txt}(46,194 real changes made)

{com}. replace NonCNRyrs2t1=0
{txt}(46,194 real changes made)

{com}. replace NonCNRyrs3t1=0
{txt}(46,194 real changes made)

{com}. replace CNRyrst1=0
{txt}(1,502 real changes made)

{com}. replace CNRyrs2t1=0
{txt}(1,502 real changes made)

{com}. replace CNRyrs3t1=0
{txt}(1,502 real changes made)

{com}. predict p1
{txt}(option {bf:pr} assumed; Pr(civicnrt1))

{com}. predict ihat, xb
{txt}
{com}. predict error, stdp
{txt}
{com}. generate lb = ihat - invnormal(0.975)*error
{txt}
{com}. generate ub = ihat + invnormal(0.975)*error
{txt}
{com}. generate plb1 = normal(lb)
{txt}
{com}. generate pub1 = normal(ub)
{txt}
{com}. drop ihat error lb ub
{txt}
{com}. 
. replace smlmegip=1
{txt}(0 real changes made)

{com}. replace highdemoc=0
{txt}(47,606 real changes made)

{com}. predict p2
{txt}(option {bf:pr} assumed; Pr(civicnrt1))

{com}. predict ihat, xb
{txt}
{com}. predict error, stdp
{txt}
{com}. generate lb = ihat - invnormal(0.975)*error
{txt}
{com}. generate ub = ihat + invnormal(0.975)*error
{txt}
{com}. generate plb2 = normal(lb)
{txt}
{com}. generate pub2 = normal(ub)
{txt}
{com}. drop ihat error lb ub
{txt}
{com}. 
. replace smlmegip=0
{txt}(47,606 real changes made)

{com}. replace highdemoc=1
{txt}(47,606 real changes made)

{com}. predict p3
{txt}(option {bf:pr} assumed; Pr(civicnrt1))

{com}. predict ihat, xb
{txt}
{com}. predict error, stdp
{txt}
{com}. generate lb = ihat - invnormal(0.975)*error
{txt}
{com}. generate ub = ihat + invnormal(0.975)*error
{txt}
{com}. generate plb3 = normal(lb)
{txt}
{com}. generate pub3 = normal(ub)
{txt}
{com}. drop ihat error lb ub
{txt}
{com}. 
. replace smlmegip=0
{txt}(0 real changes made)

{com}. replace highdemoc=0
{txt}(47,606 real changes made)

{com}. predict p4
{txt}(option {bf:pr} assumed; Pr(civicnrt1))

{com}. predict ihat, xb
{txt}
{com}. predict error, stdp
{txt}
{com}. generate lb = ihat - invnormal(0.975)*error
{txt}
{com}. generate ub = ihat + invnormal(0.975)*error
{txt}
{com}. generate plb4 = normal(lb)
{txt}
{com}. generate pub4 = normal(ub)
{txt}
{com}. drop ihat error lb ub
{txt}
{com}. 
. duplicates drop p1 plb1 pub1 p2 plb2 pub2 p3 plb3 pub3 p4 plb4 pub4, force

{p 0 4}{txt}Duplicates in terms of {res} p1 plb1 pub1 p2 plb2 pub2 p3 plb3 pub3 p4 plb4 pub4{p_end}

{txt}(47,605 observations deleted)

{com}. keep p1 plb1 pub1 p2 plb2 pub2 p3 plb3 pub3 p4 plb4 pub4
{txt}
{com}. stack p1 plb1 pub1 p2 plb2 pub2 p3 plb3 pub3 p4 plb4 pub4 , group(4) clear
{txt}
{com}. rename _stack type
{res}{txt}
{com}. label define type 1 "S Size, High Democ" 2 "S Size, Low Democ" 3 "L Size, High Democ" 4 "L Size, Low Democ"
{txt}
{com}. label values type type
{txt}
{com}. eclplot p1 plb1 pub1 type, scale(1.25) aspectratio(.65) ytitle("") xtitle("") xsize(4) xlabel(1 2 3 4)
{res}{txt}
{com}. graph save Graph figure1.gph, replace
{txt}(note: file figure1.gph not found)
{res}{txt}(file figure1.gph saved)

{com}. 
. clear
{txt}
{com}. log close
      {txt}name:  {res}<unnamed>
       {txt}log:  {res}C:\Stata\log_figure1.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res}14 Jul 2018, 19:50:08
{txt}{.-}
{smcl}
{txt}{sf}{ul off}